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Abstract 

Compact U(l) lattice gauge theory in four dimensions is studied by means of an efficient 
algorithm which exploits the duality transformation properties of the model. We focus our 
attention onto the confining regime, considering the interquark potential and force, and the 
electric field induced by two infinitely heavy sources. We consider both the zero and finite 
temperature setting, and compare the theoretical predictions derived from the effective string 
model and the dual superconductor scenario to the numerical results. 


1 Introduction 


Confinement of the color degrees of freedom in chromo-singlet hadronic states is one of the most 
intriguing phenomena taking place in the low energy regime of quantum gauge theories. The 
problem cannot be studied perturbatively, but a number of different tools have been used, in 
order to find a theoretical description of confinement; in the years, several different objects have 
been proposed as responsible for the mechanism driving this phenomenon — see [1,2] for a review 
— and, although a complete analytical solution is still missing, considerable progress has been 
achieved towards the theoretical interpretation of this problem [3]. 

Regardless of the non-universal details in their microscopic dynamics, different gauge theories 
share common properties in the infrared regime of the confined phase, where some effective degrees 
of freedom may account for the observed phenomena. 

In particular, according to the effective string scenario [4-6], the long distance properties of 
the potential of a confined charge-anticharge pair (QQ) in the heavy quark limit are expected 
to be described in terms of a fluctuating string; as it will be discussed in the following, some 
predictions of this picture have been successfully confirmed by numerical tests, but other aspects 
are still to be understood. Also, lattice results seem to indicate that the fine details could be 
different for different confining gauge theories. 

At present, the details about a possible effective theory to provide a description of gauge 
models in the IR regime are not completely known yet, and a closer look at the behaviour of 
different gauge models may be useful to get some insight into the confinement dynamics. 

In this paper, we focus our attention onto compact U(l) lattice gauge model in D = 4: an 
abelian theory which is known to possess a confined phase, analogous to non-abelian models. We 
restrict our attention to the infinite quark mass limit, considering a system with a static QQ 
pair of external charges, and work out numerical results to be compared with the theoretical 
expectations for the interquark potential V(r) and for the force F(r). We study both the zero 
and finite temperature cases. For this model, the values of the transition critical coupling at 
different finite temperatures have been published in a recent study [7], and for these values we 
focus our attention onto different ratios T/To in order to point out the quantitative effects that 
show up at finite temperature. 

Then, we study the electric field induced by the QQ pair, measuring the flux-tube profile in 
the middle plane between the two charges: according to the “dual superconductor” hypothesis, 
in the confined phase the flux lines joining the sources are expected to get “squeezed” in a tube. 
Similar studies about this model can be found in the literature — see, for instance, [8-10]. 

The numerical results are obtained by means of a Monte Carlo algorithm which exploits an 
exact, analytical mapping of U(l) gauge theory to a dual model; this procedure allows to get 
high precision results. 

The structure of the paper goes as follows: in section El the effective string scenario for 
confinement is described; next, in section El the basic properties of compact U(l) lattice gauge 
theory are reviewed. Section |11 presents the duality transformation mapping the U(l) theory 
to a model with integer variables, while in section 0 a few details about the algorithm and the 
input parameters are listed; then we discuss the way to measure the different observables, and 
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present numerical results: the latter are compared with theoretical predictions, and their physical 
meaning is discussed. Finally, in section El we summarise the main results of the work, stressing 
some concluding remarks. 

An account on preliminary results of this work has been presented at the XXII International 
Symposium on Lattice Field Theory (Lattice 2004), Batavia, IL, USA, 21-26 June 2004 [11]. 


2 Effective string theory 


The ground state interquark potential V(r) of a heavy QQ pair in a pure gauge theory (without 
dynamical matter fields) can be expressed in terms of the two-point correlation function G(r) 
among Polyakov lines: 

V(r) = -^ log G(r) = -^ log(pl(r)P(0)) (1) 

where r is the interquark distance, while L is the system size in the (compactified) “time-like” 
direction; for systems at zero temperature, the L —> oo limit is taken. 

In the confining regime, the asymptotic behaviour of V (r) at large distances is a linear rise, 
a fact compatible with the idea that the chromoelectric flux lines get squeezed into a “tube” 
between the two charges, whose energy is expected to be (roughly) proportional to its length. It 
was therefore conjectured that the long distance properties characterising the heavy QQ system 
in the confinement phase can be modelled in terms of a vibrating, string-like object joining the 
two sources; this approach, which nicely fits the intuitive picture of a chromoelectric flux tube 
squeezed into a thin Abrikosov-Nielsen-Olesen vortex [12,13], was quantitatively developed by 
Liischer, Munster, Symanzik and Weisz [4-6] more than two decades ago. 

Obviously, in a realistic model, matter fields would be dynamically created out of the vacuum 
and induce the “string breaking” phenomenon, leading to the flattening of the interquark potential 
at large distances; however, here we restrict our attention to the case of a pure gauge theory, with 
no dynamical matter fields. 

According to the effective string scenario, the IR behaviour of V(r) should be predicted from 
the form of the action governing the dynamics of the string-like flux tube between the static 
sources. In general, the contributions beyond the linearly rising term in V(r) are expected 
to depend on the functional form of the effective string action, but the coefficient of the 1/r 
contribution (which is known as the “Liischer term”, the leading string correction in the IR 
regime) can be worked out under the very general assumption that, at large distances, the effective 
model reduces to a free vibrating string. A first quantisation approach shows that the leading 
contributions to V(r) in the string picture read [4]: 


V(r) 


ar + /j, 
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where D is the number of spacetime dimensions, a is the string tension, and /r is a constant. 
The Liischer term in eq. m must not be confused with a Coulomb term originating from the 
one-gluon exchange process (which is relevant at short distances); on the contrary, the Liischer 
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term is a Casimir effect due to the finiteness of the interquark distance r, and it is a consequence 
of the hypothesis that, at leading order, the transverse string fluctuations are described by a 
CFT of {D — 2) free, massless boson fields. Actually, it was shown that —is a universal 
coefficient that one expects in any model based on a string theory [6,14], 

Whether the effective string model provides a satisfactory picture or not, is still a question 
under debate. On one hand, the interquark potential V{r) (and related quantities) have been 
investigated in a large number of lattice studies for different gauge theories since the Eighties — 
see [15,16] for references — with results confirming the bosonic string prediction: in particular, 
the Liischer term appears to be a universal feature of the IR regime of confined gauge theories. 
On the other hand, data from recent, high precision Monte Carlo simulations [10,17-26] have 
shown that higher order corrections to V (r) might be non-universal, and — in general — do not 
trivially fit the bosonic string predictions. 

A related aspect, which is still under debate, concerns the excitation spectrum; according 
to the bosonic string picture, it is expected to be described by a tower of harmonic oscillator 
excitations with energies: 

E = Eo + —n , n G N (3) 

r 

where degeneracies are expected for re > 1. However, while the Liischer term in the ground state 
potential appears to be confirmed by numerical evidence down to surprisingly short distances, 
the status of lattice simulation results for the excited states seems to be much more difficult to 
fit into this simple harmonic oscillator pattern: in some studies [27-30] it was found that, at 
distance scales around 0.5 to 1.0 fm (where the Liischer term in the ground states is observed) 
the excitation spectrum of the static interquark potential does not agree with eq. ©• A similar 
mismatching between effective string predictions [31] and numerical results was also observed for 
the torelon states [32]. 

Another prediction of the effective string models concerns the behaviour of the “flux tube 
width” as a function of the distance between the sources [5]: in the confinement phase, the gluon 
field in the ground state of a heavy QQ pair is essentially expected to be non-vanishing inside a 
tube-like region between the charges, and to be exponentially decreasing with the distance from 
the interquark axis. A measure of the (square) typical transverse width of the flux tube can be 
expressed as: 

2 f d^x±xj_£(x) 

^ f d‘^x±£(x) 

evaluated in the mid-plane between the charges, where x_l is the transverse coordinate, while £(x) 
is the chromoelectric field energy density of the QQ ground state with respect to the vacuum: 

£(x) (X (QQItrE^(x)IQQ} - {QQ\QQ){^\tiE‘^{x)\Q) (5) 

In the effective string scenario, the flux tube width measures the amplitude of string vibrations, 
and, in the Gaussian approximation, it is expected to increase logarithmically as a function of 
the interquark distance: 

= ( 6 ) 
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(where A is a length scale, and Wq is a constant); this behaviour was indeed confirmed by numerical 
simulations of Z 2 gauge model in Z) = 3 [33] and U(l) gauge model in Z) = 4 [9]. 

Focusing the attention onto the ground state interquark potential, the prediction for V{r) 
of a QQ pair eq. 0 can be refined as it follows. Assuming that the spectrum (including 
multiplicities) of the pure gauge model is described by the effective string, and neglecting excited 
state decays, the partition function for this sector of the theory is given by the Polyakov loop 
two-point correlation function G{r). The latter is written as a string partition function: 

G{r) = {P\r)P{0)) = j [Vh]e-^-« (7) 


where S^s is the effective action for the world sheet spanned by the string, and functional inte¬ 
gration is done over world sheet configurations having the Polyakov lines as a fixed boundary. 
At leading order, one can assume that the degrees of freedom associated with transverse string 
vibrations are free, massless bosons^, in agreement with the physical content of eq. ©• 

To get more detailed information about the subleading terms that do not appear in eq. 0, 
one should know the higher-order contributions to the effective string action, which introduce 
string interactions. 

Choosing to be proportional to the area spanned by the string world sheet one gets the 
Nambu-Goto action [35-37]: 

S’eff = a • y d^^v^deF^ (8) 

This string action suffers from inconsistency problems at a quantum level: rotational invariance 
is broken once one considers the problem from the quantum point of view, unless the number of 
spacetime dimensions is ZA = 26 [38]. 

In the physical gauge‘s the integrand appearing in eq. (0 can be written as: 

[1 -F {doh)‘^ + (dih)'^ + {doh x (9i/i)^] ^ (9) 


(h being the transverse displacement of the world sheet surface). 

Although the resulting theory is non-renormalizable, expanding eq. © in a perturbative 
series, the (regularised) expression for (Pf(r)P(O)) at leading order reads: 


^—arL—fiL 

where Dedekind’s r] function: 


(pt(^)P(O)) 
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T/(r) = ~ ’ with: q = (11) 

_ T?.= l 

^This hypothesis is justified in the so-called “rough” regime [34] of the theory. 

^Since eq. 0 is invariant with respect to reparametrization and Weyl transformations, a gauge choice is needed 
to proceed to quantum treatment. In the physical gauge, two degrees of freedom are associated with the “natural” 
coordinates, and the remaining ones describe transverse displacements. We neglect the problems related to Weyl 
anomaly, since the latter vanishes at large distances. 
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encodes the contribution of massless fluctuations. Notice that — as it is expected — eq. m 
reproduces the Liischer term contribution to the interquark potential in the “low temperature 
regime” ^ > 1- 

The exact spectrum of the Nambu-Goto string can be obtained via formal canonical quanti¬ 
sation [39] (see also [40]); the energy levels for a string with fixed ends read; 


En{r) 
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The string partition function can be written in terms of these energy levels as: 

CXD 

Z = (13) 

n=0 

with some integer weights Wn, which account for state multiplicities. The simplest example is the 
D = 3 case, where the Wn coefficients are the partitions of n: 


Z = y]P(n)e“-®"W^ (14) 

n=0 


which has been recently compared with numerical results of Z 2 , SU(2) and SU(3) gauge theories 
[26]. 

The picture of an effective string can be generalised in different ways: for instance, in [41], 
Polchinski and Strominger proposed a non-polynomial action: 
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(15) 


(where we followed the notation of the original paper); here, a denotes the typical length scale 
related to the string tension. In H = 4, eq. (USD describes an interacting, Poincare-invariant 
CFT with four bosons and central charge c = 26; the second term appearing in the integrand 
on the right hand side of eq. m is related to the determinants discussed by Polyakov in [42], 
although here they are meant to be built out of the induced metric instead of the intrinsic one. 
More precisely, the idea underlying the construction of the effective action given by eq. m is to 
allow a generic (but conformally invariant) world-sheet QFT, and to absorb the Jacobian into the 
action coefficients; the latter can be fixed in such a way that the conformal anomaly is cancelled 
in any dimension D. The non-polynomial terms in eq. m are not pathological, provided one 
considers an expansion about the long-string vacuum^. It has been recently pointed out [43] that 
Poincare invariance also constrains the 0(r~^) term of the string spectrum to be exactly the same 
that shows up in the Nambu-Goto model. 

®When the string length is much larger than a, the infinite number of terms originating from the expansion of 
eq. 1151 are irrelevant, as they are suppressed by powers of a. 
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Eq. ((SI) and eq. m can be considered as particnlar examples for a candidate effective string 
action; more generally, the effective string action is expected to be written as a series like [18,44]: 

= crrL + + So + Si + S 2 + ..., with: Sq = ^ J (f^{dahdah) (16) 

S'o, the leading term (beyond the “classical” contributions), describes a conformal model, while 
the subsequent Sn terms, with n > 1, have increasing dimension, and introduce string self¬ 
interactions. 

In particular, 5i is associated with a “boundary contribution” [18]: 

^1 = 1 / [{dih)l^o + idih)l=r] (17) 

b being a coupling with length dimension. A perturbative calculation [23,44] shows that the effect 
of Si amounts to a shift r —> r — 6 in the contribution induced by the So term^. The role possibly 
played by this “boundary contribution” has been studied in detail in some recent works: on the 
numerical side, comparison of simulation results for different gauge models D = 3 and D = A 
seems to suggest that the 5i contribution is not universal [10, 23, 24, 30] — as opposite to the 
Liischer term induced by So- On the theoretical side, it has been recently pointed out [44] that 
the O (r“^) contribution induced by Si has to be ruled out, if open-closed string duality holds. 

The effective string scenario described here refers to a “meson-like” QQ system, but it can be 
generalised to the case of “baryonic” structures of three quarks as well [45]; however, here we do 
not address this issue, restricting our attention to predictions for the QQ system that we studied 
in lattice simulations. 

3 Compact U(l) lattice gauge theory 

We study compact U(l) gauge theory in four spacetime dimensions, on a hypercubic, isotropic, 
Euclidean lattice A with spacing a; periodic boundary conditions are imposed along all directions; 
here and in the following, denotes the number of lattice sites along the fi direction. 

This simple abelian theory has some interesting properties in common with non-abelian models 
— remarkably, the existence of a confined phase — and, historically, it has been the object 
of many numerical studies since the early days of lattice gauge theory. The analogy with its 
three-dimensional counterpart, where charge confinement mechanism is driven by condensation 
of monopole plasma [46,47], is supported by numerical evidence that the monopole density is 
non-vanishing in the confined phase (at strong coupling), while it goes to zero in the deconfined 
regime. 

The fundamental degrees of freedom of the model are U^{x) = exp [i9^{x)] e U(l) variables 
defined on the oriented links of the lattice; their dynamics is governed by the standard Wilson 
action: 

S = (3 [1 — cos {dO)] (18) 

_ pi- 

■^In particular, this holds up to and including 0{b^) terms. 
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where the summation ranges over the lattice plaquettes, d denotes the discretised exterior deriva¬ 
tive, and (3 = 

At zero temperature, the theory has two different phases: for f3 values below a critical value 
Pc = 1.0111331(21) [48], the system is in the confined phase, whereas for > Pc it lies in a 
Coulomb-like phase; the transition at P = Pc is a (weak) first order one. In the present work, we 
consider values of P in the confined phase, close to the deconfinement transition®. 

As it concerns the finite temperature case, in [7] the phase diagram of the theory was studied 
for different values of the number of Euclidean time slices, using the helicity modulus — see 
also [49] — as a tool to investigate the phase transitions; the numerical findings show that in the 
{P, ■^) plane a phase boundary separates the confining phase from a temporal Coulomb phase 
for any finite value of A4. Moreover, the critical values of the coupling parameter Pc at which 
the system undergoes a finite temperature phase transition, and the corresponding values of the 
latent heat were measured; the latter turns out to be a quantity which is finite for large values of 
Ni (a signal of the first-order nature of the phase transition at T = 0), but it quickly decreases 
for growing temperatures. In particular, for A4 < 4, the latent heat appears to vanish, opening 
the possibility of a second order phase transition at high temperatures. 

The interquark potential for a QQ system can be expressed in terms of the expectation value 
of the correlator among Polyakov lines P{x) as: 

V{r) = \og{P\r, 0, 0)P(0,0, 0)) (19) 

We also study the behaviour of the longitudinal component of the electric field E in the \QQ) 
state; E is expressed in terms of the field strength F, which can be defined as; 

F = -^^sin((i0) (20) 

This definition is consistent with the electric Gauss law. 

4 Duality transformation 

Compact U(l) lattice gauge theory can be analytically mapped to a dual formulation via a group 
Fourier transform: the partition function and expectation values of the original theory can be 
rewritten in terms of variables of the dual model. 

This technique is inspired by spin duality in statistical systems [50,51]; in its simplest formu¬ 
lation, it can be used for an arbitrary abelian or solvable group Q [52, 53] in generic spacetime 
dimension D. This is due to the fact that, for abelian groups, the product of irreducible rep¬ 
resentations still gives an irreducible representation; for solvable groups, the reasoning can be 
generalised considering a well-suited mapping of group elements down to an abelian group — see 

®Since the deconfining transition is first order, the correlation length scale ^ = (j“ 2 is not divergent as the critical 
point is approached, thus, strictly speaking, one cannot expect to exactly recover the continuum limit; however, 
the limit value of the string tension for /? —> /Jc appears to be quite small (approximately of order 
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also [54]. On the other hand, a generalisation to non-abelian lattice gauge theory was discussed 
in [55], but in that case the resulting spin foam model involves non-trivial constraints. 

Among the main features of this duality mapping, it is worth stressing its intrinsically non- 
perturbative nature, and the fact that it is not a configuration-to-configuration mapping. 

Besides being theoretically interesting, this mapping also proves to be useful in numerical 
studies, since evaluation of physical observables from measurements of related quantities in the 
dual model often offers practical advantages. Improvement may occur due to various reasons: the 
dual model may be defined in terms of variables which are “simpler” (or “simpler to handle” for 
a computer) than the original ones; the dynamics of the dual model itself may be “simpler” than 
in the original theory; some observables may be evaluated with a higher precision from direct 
measurements in the dual model. Indeed, the duality technique has been successfully used in 
many numerical studies of Z 2 lattice gauge theory [19,23,30,56-58], and it was already tested for 
this model in [9], where the electric field, magnetic current, and flux tube width were investigated; 
on the other hand, in [59-61] it was applied to the Villain formulation of the theory. 

As a matter of fact, the duality transformation allows one to trade the continuous variables 
of U(l) model for integer-valued variables, but the major improvement can be achieved in the 
study of the interquark force as a function of distance. In the original model, the interquark 
force is expressed in terms of ratios between expectation values of Polyakov loop correlators, 
that are exponentially decreasing with the distance r; these quantities are thus affected by a fast 
signal-to-noise ratio decay at large distances, and — in general — error-reduction techniques are 
needed to get high precision results. In the dual formulation, instead, the problem is completely 
overcome, since a lattice estimate of the interquark force for the U(l) gauge model — which, in 
se, is finite even for very large interquark distances — can be obtained from direct measurements 
in the dual system. In this respect, this duality-inspired method can be considered as a possible 
alternative to other popular error-reduction numerical techniques, like the multi-level algorithm 
proposed by Liischer and Weisz [62], which has been used to achieve exponential error reduction 
in several different gauge models [18,20,24,63,64], including compact U(l) lattice gauge theory 
in D = 4 [10,65]. 

The duality transformation for compact U(l) lattice gauge theory in 11 = 4 can be summarised 
as it follows (see also [9]): the plaquette terms appearing in the partition function: 

Z = [ Y]. I [1 “ (d.0)] > (21) 

x,n I pi. J 

are expanded in a group Fourier series: 

^ = / Y{def,{x)Y{'^exp[i{k,de)]e~^Iik\iP) ( 22 ) 

x,n pi. fcez 

/|fc|(z) is a modified Bessel function of the first kind of order \k\, and k is an integer-valued 2-form 
on the lattice, which is associated with plaquettes. Rewriting the product (A:, dO) as: {5k, 9), and 



integrating over the 9 link variables introduces a constraint on the /c’s: 


Z = (27r)4^nj;e-%|(/3) 

pi- /cGZ 


Sk=0 


(23) 


(4A^ being the total number of links). Physically, k is associated to the field strength, and the 
constraint 6k = 0 enforces the vacuum Maxwell equation: 


dF = Je = 0 

The constraint 6k = 0 is solved introducing the integer-valued 3-form I such that: 


(24) 


k = 6l 


(25) 


and the partition function, written as a sum over configurations (to be denoted as X){i}) 
“elementary cube” variables I taking values in Z: 


Z = (27r)4^^n 

{/} pi- 


(26) 


can be reformulated in a dual lattice description, where the I 3-form defined in the direct lattice 
corresponds to a *l 1-form living on the oriented links of the dual lattice *A: 


Z = (27r)4^^n[e-%*,|(/3) 

{* 1 } pi- 


(27) 


(in this case, the product ranges over the plaquettes of the dual lattice). Eq. describes a 
model whose fundamental degrees of freedom are *l link variables interacting through a plaquette- 
like term: 

Z = ^ exp [-5dua,i] (28) 

{*1} 

which is a sum over all the configurations of integer-valued *l variables defined on the dual lattice 
links, with: 

*8^dual E [P - log I\dn\ (/3)] + const (29) 

pi. 

The expression for S'duai is thus given in terms of a sum over the plaquettes of the dual lattice, 
and it shows that the dual model is again a “gauge model”, because S'duai only depends on d*l 
(thus it is left invariant under integer-valued local gauge transformations, defined on the *A sites, 
that shift the incoming */’s by the same amount) and, since IniP) is a decreasing function of the 
integer n index, for any fixed, positive value of P, this interaction has “ferromagnetic nature”, 
in the sense that the configurations with the largest statistical weight are the ones making \d*l\ 
minimal. 
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An analogous expression can be derived when external charges are included; in that case, the 
starting point is the QQ system partition function; 

r ( 1 ^4-1 N 4-1 

ZqQ= U dW E [1 - (de)] n n (30) 

x,fi t pl- J T=0 

which is related to the Polyakov loop two-point correlation function: 

(pt(,)p(o)) = ^ (31) 

In complete analogy with Z, eq. (|30|) can be rewritten in the dual formulation. In this case, 
the constraint 5k = 0 modifies to 5k = ±1 for links belonging to the (anti)charge Polyakov loop, 
and this can be enforced introducing a further 2-form n which describes a world sheet having the 
Polyakov lines as its boundary: 

k = 51 + n (32) 

and which is related to the presence of non-vanishing external charges: 

e5n = J (33) 

The precise (local) shape of the stack of non-vanishing re’s is irrelevant, and the simplest choice 
corresponds to the minimal surface tiling the sheet between the Polyakov lines — see also fig. ^ 
showing the *re stack of defects as shaded plaquettes in the dual lattice. However, in principle 
one should keep into account all of the different global (topological) structures for the stack of 
defects, and include possible Z-twists in the boundary conditions, which create interfaces. The 
cartoon in fig. |21 shows a two-dimensional sketch of two different topologies for the *re stack: on 
the right hand side, an interface is created. The larger energy cost associated to this topology 
can be compensated by the gain in entropy, due to the fact that the interface can be created 
anywhere. However, non-trivial topologies turn out to be statistically suppressed in the confining 
regime, for large enough lattice cross-sections, and in our numerical simulations we neglected 
their effect. 

The expression for Zqq in the dual lattice reads: 

Zqq = (27r)'^^ E I “ E I\dH+*n\ (/?)] } (34) 

{*1} I pi. ) 

the only difference with respect to eq. (EHl) being the presence of the *re defects which account 
for the static charges. 

Finally, to study the field strength flux through a given plaquette p in the QQ system, one 
further needs to include the y^sin {dO)^ operator in the path integral; the expression for {Fp)qq 
can be worked out using the following property of modified Bessel functions of the first kind: 

2?7- 

—In{z) = In-l{z) - In+l{z) (35) 
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Figure 1: The QQ pair, described by Polyakov loops in the original model (solid lines), is mapped 
to a stack of topological defects {*n = +1) localised onto the shaded plaquettes in the dual lattice. 
In particular, the shaded plaquettes are dual to an (arbitrary) surface in the original lattice, having 
the Polyakov lines as its oriented boundary. 
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Figure 2: Schematic, bidimensional, view of two different topological structures for the *n stack 
of defects (solid black lines) in a system with torus-like geometry. On the left side, the defects 
are located along the shortest line joining the (point-like) sources: this is the analogue of the 
scheme depicted in hg. ^ On the right side, the line of defects joining the sources winds around 
the system, and a further interface is generated. The latter, by virtue of translational invariance, 
accounts for an “entropy gain”, which can compensate for the larger energy cost with respect to 
the pattern shown on the left side of the figure. 
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and the result for {Fp)qq reads: 


i^p)QQ = 7 — ( 2 vr)'^^ ^ I “ I\dH+*n\ (/3)] > (36) 

I Pi. J 

5 Numerical simulations: settings and results 

The algorithm performs Monte Carlo simulations of the dual lattice model, updating the *l link 
variables; the *n variables that account for the presence of external charges are kept fixed. We 
chose the simplest *n configuration, like in figure ^ checking that, for the given values of the 
parameters, the effect of different topologies could be neglected within the precision of our results. 

In the runs, we considered lattice sizes up to 19 x 31^; the /? coupling — the parameter of 
the original theory — took values from 0.96 to 1.0107, in the confined regime and close to the 
critical value jdc- The configuration update was done via standard Metropolis method; this choice 
highlights the advantages arising from the use of the duality transformation. 

The autocorrelation time was evaluated using the automatic windowing procedure due to 
Madras and Sokal [66]; then we chose to collect sets of non-autocorrelated data only. The errorbars 
from data samples have been obtained via the jackknife method. 

Detailed results are presented in the following subsections. 

5.1 Interquark potential and force at zero temperature 

The interquark potential V{r) and force F{r) can be obtained from the expectation value of the 
two-point Polyakov loop correlators; we run zero temperature simulations on a hypercubic lattice 
of size 16^, considering different (on-axis) interquark distances. The values of the j3 parameter 
were chosen in the range between 0.975 and 1.010: a window in the confined phase, close to 
the deconfinement critical point jdc] correspondingly, the value of the string tension a ranges, 
approximately, from 0.3a“^ down to 0.06a“^, so that the product crL^ is large enough to damp 
the statistical weight of configurations with an interface (see discussion in sect. 0 and fig. EJ. 
This choice for the parameters also allows a direct comparison with the results obtained with a 
different algorithm in [10]. 

In particular, it is interesting to note that our algorithm allows to measure easily the ratio 
between Polyakov loop correlators at distances r -|- 1 and r: in the dual lattice, such a ratio can 
be obtained in terms of one column of plaquettes, in the background of a L x r stack of defects 
describing the QQ pair at distance r. This enables us to overcome the problems due to the fact 
that G{r) decreases exponentially fast with r: although both G{r + 1) and G{r) are quantities 
which become very small already at distances of a few lattice spacings, their ratio goes to a finite 
limit for large r, and it is expected to develop (effective string) corrections for finite values of 
r. In order to achieve a further improvement in the result precision, we actually split the ratio 
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appearing on the right hand side of eq. m into a product of simpler ratios: 


G{r) 


^QQ 

Z 


n 


^{Lxr)+i+l 

^{Lxr)+i 


( 37 ) 


where the generic term ^(i,xr)+i appearing on the right hand side of eq. (1371) describes the partition 
function of a system where a column of i defects is attached to the original L x r stack, so that 
^{Lxr)+i .Z’(^xr)+i+i differ for the presence of a *n defect on one plaquette. Keeping 
eq. ( 1511 ) in mind, it is then easy to see that the factors appearing on the right hand side of eq. (| 57 |) 
are nothing but the expectation values of the local quantity: 


|d*/|(/3) 


evaluated in a system with a (L x r)+i stack of *n defects. This is the so-called “snake algorithm” 
idea, which was introduced in [67,68], and that has also been successfully used in [19]. Each of 
the factors appearing on the right hand side of eq. (| 57 |) . which are obtained from separate, 
independent simulations, are affected by small variance. 

The first task is to check the presence of the leading-order (LO in the following) correction 
with respect to the asymptotical linear rise of V{r) which is expected at large distances. In the 
effective string scenario, in H = 4 and at zero temperature, eq. GDI) gives: 


V{r) = ar + fj,- + O (r 


(39) 


Eq. (irni) can be fitted to the numerical results for G{r). 

Results of data analysis for the different values of /3 that were investigated are reported in 
tab. n since eq. m predicts the behaviour at large distances, we progressively discarded the 
G{r) data corresponding to the smallest r values. 

The results obtained for the string tension a are plotted as a function of f3 in fig. |5J 

In tab. [21 we also report the results that we obtained for Sommer’s scale ro [69], which is 
expected to provide a natural scale to express the numerical results; as an example, in fig. [^ the 
plot for the interquark potential is shown: the good scaling properties of the data obtained from 
different values of f5 are manifest. 

As it concerns the interquark force F, evaluated at the mid-point between r and r-|-l, note that 
it can be approximately estimated in terms of the Polyakov loop correlator ratio G{r + 1)/G(r) 
as: 


F 






G(r + 1) 
G{r) 


(40) 


As we already remarked, the results for G{r + 1)/G{r) that can be obtained by means of our 
algorithm at larger and larger distances r are not affected by the typical (exponential) signal-to- 
noise ratio decay, which is observed for Monte Carlo simulations in the direct model; this is a 
major advantage that comes from simulations in the dual setting, as it was also pointed out in 
the study of the Z 2 gauge theory in D = 3 in [19]. 
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Figure 3: Results for the string tension a (in lattice units), as a function of (5. The dashed line 
corresponds to the critical value (3c = 1.0111331(21) [48]. 
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Figure 4; The interquark potential as a function of the distance, in units of Sommer’s scale tq. 
Different symbols refer to different values of (3, in the range from 0.975 to 1.010. 
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/3 


cra^ 

/iO 

number of d.o.f. 

reduced 

0.975 

1 

0.31322(95) 

0.4996(33) 

6 

0.8 


2 

0.3145(13) 

0.4929(59) 

5 

0.7 


3 

0.31208(50) 

0.5068(25) 

4 

0.6 

0.980 

1 

0.28833(79) 

0.5114(29) 

6 

1.3 


2 

0.2889(12) 

0.5082(55) 

5 

0.9 


3 

0.2869(10) 

0.5196(53) 

4 

0.8 

0.985 

1 

0.2616(15) 

0.5189(51) 

6 

1.7 


2 

0.2608(24) 

0.523(10) 

5 

1.9 


3 

0.25671(97) 

0.5453(50) 

4 

0.7 

0.990 

1 

0.2277(12) 

0.5273(42) 

6 

0.6 


2 

0.2287(18) 

0.5226(83) 

5 

0.7 


3 

0.2256(14) 

0.5398(74) 

4 

0.3 

0.995 

1 

0.20085(62) 

0.5447(22) 

6 

0.4 


2 

0.202128(60) 

0.5379(27) 

5 

0.2 


3 

0.20184(91) 

0.5396(47) 

4 

0.2 

1.000 

1 

0.16357(66) 

0.5550(24) 

6 

0.4 


2 

0.1634(10) 

0.5561(46) 

5 

0.4 


3 

0.16349(88) 

0.5570(49) 

4 

0.3 

1.005 

1 

0.12252(73) 

0.5718(27) 

6 

0.5 


2 

0.1231(11) 

0.5686(49) 

5 

0.5 


3 

0.12302(90) 

0.5693(33) 

4 

0.8 

1.010 

1 

0.06149(87) 

0.5798(32) 

6 

0.6 


2 

0.0627(11) 

0.5699(51) 

5 

0.7 


3 

0.06291(67) 

0.5734(29) 

4 

0.6 


Table 1: Fit of eq. HIOI) to the simulation results for the Polyakov loop correlator G{r). Data 
corresponding to distances below are discarded. 


Tab. 131 displays our results for G{r + 1)/G(r) at different values of r and j3. Notice that, for 
each set, the uncertainty associated to the data at large values of r appears to be fairly constant, 
or growing slowly; since the number of measurements to obtain these results was the same for all 
distances in each set, this improvement in precision with respect to the standard approaches can 
be fully ascribed to the properties of the algorithm. 

According to eq. m, the LO, large-distance behaviour which is expected for the force reads: 

+ ( 41 ) 

Expressing our results in terms of Sommer’s scale ro, it is possible to make a comparison of the 
whole set of data obtained from the simulations at values of (3 ranging from 0.975 to 1.005 with 
the LO prediction obtained from the effective string picture. The corresponding plot is shown 


16 


















3 

2.5 

2 

ti, 


1 

0.5 

0 

0 0.5 1 1.5 2 2.5 3 3.5 

r/rg 

Figure 5; The approximate interquark force F — obtained from finite (one-lattice-spacing) dif¬ 
ferences of the potential — as a function of the (mid-point) distance. All data, obtained from 
simulations at different values of f3 in the range from 0.975 to 1.005, are rescaled in terms of 
Sommer’s scale ro (which was obtained from separate sets of simulations). The dashed line is the 
theoretical expectation for the interquark force, as it is predicted by the effective string picture 
at LO. 
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p 

ro/a 

0.975 

2.1091(17) 

0.980 

2.1997(38) 

0.985 

2.3254(44) 

0.990 

2.4806(77) 

0.995 

2.6225(60) 

1.000 

2.9147(89) 

1.005 

3.358(15) 

1.010 

4.751(34) 


Table 2; Results for Sommer’s scale tq at different values of f3. 

in fig. El Note that, at this level, the dashed curve does not depend on any free parameter, yet 
it fits quite well with the numerical data, despite the approximations that we discussed. The 
matching between LO effective string prediction and numerical results appears to hold down to 
short distances — a fact in agreement with similar conclusions that have been obtained in other 
studies of different gauge theories [18,70]. 

5.2 Finite temperature setting 

The analysis of subsection t1.ll can be extended to the finite temperature case. There are many 
physical motivations for addressing the study of gauge theories at finite temperature, and, in 
particular, it is remarkably interesting to see what happens as one considers temperatures which 
are closer and closer to the finite temperature deconfinement transition: in this subsection, we 
present the results for the Polyakov loop correlator ratios at finite temperatures, up to about 3/4 
of the critical deconfinement temperature. 

The interquark potential V{r) can still be expressed by eq. (HUll . and the LO bosonic string 
prediction for the Polyakov loop two-point correlation function is given by eq. (cni), but an 
important difference with respect to the zero temperature setting is that, depending on the 
values of r and L, at finite temperatures it is possible to distinguish between two different regimes 
[17,19,24]: 

1 . ^ > 1, that we shall call the low temperature regime, 

2 . ^ < 1, to be denoted as the high temperature regime. 

The functional form of the leading term due to string effects in the ^ 1 limit of the low 

temperature regime obviously reproduces the “Liischer term” correction to V(r) appearing in the 
zero temperature case; on the other hand, in order to study the opposite ^ S> 1 limit of eq. m, 
we make use of the following property of Dedekind’s function: 

r? = V-ir r]{T) (42) 


18 






/3 

rfa 

G{r + 1)/G{r) 

P 

rfa 

G(r + 1)/G(r) 

P 

rfa 

G{r + 1)/G(r) 

0.975 

0 

1.440(80) X 10-4 

0.990 

0 

3.66(30) X 10-4 

1.005 

0 

9.71(57) X 10-4 


1 

1.002(60) X 10-3 


1 

3.93(34) X 10-3 


1 

1.90(11) X 10-2 


2 

2.70(17) X 10-3 


2 

1.069(91) X 10-2 


2 

6.64(39) X 10-2 


3 

4.49(28) X 10-3 


3 

1.74(10) X 10-2 


3 

9.62(57) X 10-2 


4 

5.61(34) X 10-3 


4 

1.99(17) X 10-2 


4 

1.012(60) X 10-1 


5 

6.13(38) X 10-3 


5 

2.57(22) X 10-2 


5 

1.361(80) X 10-1 


6 

6.16(38) X 10-3 


6 

2.63(16) X 10-2 


6 

1.317(77) X 10-1 


7 

6.15(38) X 10-3 


7 

2.63(16) X 10-2 


7 

1.522(89) X 10-1 

0.980 

0 

1.80(15) X 10-4 

0.995 

0 

4.20(24) X 10-4 

1.010 

0 

2.33(14) X 10-3 


1 

1.39(12) X 10-3 


1 

5.54(33) X 10-3 


1 

4.94(29) X 10-2 


2 

4.34(37) X 10-3 


2 

2.02(12) X 10-2 


2 

1.642(96) X 10-1 


3 

6.46(35) X 10-3 


3 

2.52(15) X 10-2 


3 

2.42(14) X 10-1 


4 

8.68(75) X 10-3 


4 

3.29(20) X 10-2 


4 

3.24(19) X 10-1 


5 

8.83(76) X 10-3 


5 

3.40(20) X 10-2 


5 

3.46(20) X 10-1 


6 

8.63(74) X 10-3 


6 

3.73(22) X 10-2 


6 

3.73(21) X 10-1 


7 

1.110(96) X 10-2 


7 

3.96(24) X 10-2 


7 

3.66(21) X 10-1 

0.985 

0 

2.53(15) X 10-4 

1.000 

0 

6.73(40) X 10-4 




1 

2.15(13) X 10-3 


1 

9.58(58) X 10-3 




2 

5.65(34) X 10-3 


2 

3.10(19) X 10-2 




3 

1.074(65) X 10-2 


3 

5.42(32) X 10-2 




4 

1.298(79) X 10-2 


4 

6.01(36) X 10-2 




5 

1.485(90) X 10-2 


5 

6.21(37) X 10-2 




6 

1.490(90) X 10-2 


6 

7.25(43) X 10-2 




7 

1.72(10) X 10-2 


7 

7.29(43) X 10-2 




Table 3: Ratios of Polyakov loop correlators at different distances and values of (5. These results 
are obtained from zero temperature simulations on 16*^ lattices. 


Accordingly, the LO string prediction for V{r) in the high temperature regime reads: 

V{r) = (^cj-^)r + ^ + ^log^ + |-^log(l-e“^’^) (43) 

n=l 

from which one can see that the string fluctuations act to renormalize the string tension, through 
the — 3^7 term. 

In the following of this subsection, we shall focus our attention on the 2r > L case: a 
regime where numerical studies for the Z 2 gauge model in D = 3 [19] have shown that the finite 
temperature effects appear to be consistent with the predictions of the Nambu-Goto effective 
string action truncated at NLO. 

In [7], the values of the critical coupling (3c for compact U(l) lattice gauge theory in = 4 
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have been obtained, for different values of the inverse temperature L. In particular, here we 
consider the L = 8 and L = G cases, for which the values of (3 at which the finite temperature 
deconfinement transition occurs are, respectively, (5c = 1.0107(1) and f5c = 1.009449(1); for 
these values of the coupling parameter, we run numerical simulations at different temperatures, 
approaching the deconfinement temperature from below. The lattice size along the remaining 
directions Ng (s = 1,2,3) is large enough, so that neglecting non-trivial topologies for the stack 
of *n defects is still a justified approximation. 

We measured the Polyakov loop correlator ratios G{r + 1)/G{r) with the same method dis¬ 
cussed in subsection 15.1 1 information about the parameters in these finite temperature simulations 
is summarised in tab. [I] The corresponding plots are shown in fig. Eland in fig. [T) 


p 

l/Tc 

l/T 

T/Tc 

Ns 

r/a range 

aa^ 

reduced 

1.009449 

6 

8 

0.75 

21 

5-9 

0.06196(91) 

0.16 



10 

0.60 

23 

6-10 

0.0665(27) 

1.73 



12 

0.50 

25 

7-11 

0.0683(25) 

1.76 



14 

0.43 

27 

8-12 

0.0662(17) 

1.02 

1.0107 

8 

11 

0.73 

23 

6-10 

0.0409(25) 

1.74 



13 

0.62 

25 

7-11 

0.0474(15) 

0.71 



16 

0.50 

29 

9-13 

0.0415(28) 

3.15 



19 

0.42 

31 

10-14 

0.0476(10) 

0.47 


Table 4: Parameter information for the finite temperature runs, together with the results obtained 
from the fit of the behaviour predicted by eq. m to the G{r + l)/G{r) ratios. 

The results displayed in tab. |1] are essentially in agreement with the behaviour predicted by 
eq. m for the two-point Polyakov loop correlation function, and the values obtained for the 
string tension have relative errors between 1% and 3% (with the exception of a unique data 
sample). 

However, data analysis also shows that the precision of these measurements at finite tem¬ 
perature does not allow to give a conclusive answer about the fine effects predicted assuming a 
string-like behaviour. 

5.3 Electric flux profile 

In this subsection we discuss the results for the electric field E in the QQ system. With no loss of 
generality, the QQ axis can be chosen to lie along the x axis of the lattice; we study the electric 
field longitudinal component Ex in the {y, z) mid-plane between the charges. 

As opposite to the Coulomb-like regime, where Ex is expected to fall off like a power of 
p = ill the confined phase the dual superconductor picture predicts an exponentially- 

driven decay which, at large enough distance from the QQ axis, reads [10]: 

Ex {p) = m%Ko {ruBp) (44) 
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G(r+1)/G(r) G(r+1)/G(r) 
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Figure 6: Fit of the theoretical expectation derived from eq. cni) to the numerical results for 
Polyakov loop correlator ratios, as a function of the interquark distance (expressed in terms of 
Sommer’s scale) at /3 = 1.009449; the corresponding critical deconfinement inverse temperature 
is 1/Tc = 6 [7]. 
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Figure 8: Electric flux through plaquettes in the mid-plane between the charges, as a function of 
the distance p from the QQ axis, at (5 = 0.96, on a 8“^ lattice. 



where niB the mass of the dual gauge boson. Eq. (jlH) is derived in the approximation of an 
infinitely long flux tube joining the sources, neglecting finite size effects and string fluctuations. 

In these simulations, the interquark distance r was kept hxed (r = 3a), and different values 
of (3 were chosen, in the range from 0.96 to 1.01; different lattice sizes 8“^ and 16^ were also 
considered, in order to check the influence of possible finite lattice size effects. Information about 
the parameters is summarised in table [SJ 


lattice size 

8^ 16^ 

P 

0.96, 0.97, 0.98, 0.99, 1.00, 1.01 

interquark distance r 

3a 

number of measurements 

20 000 


Table 5: Information about the simulation parameters in the study of the electric flux profile. 

Our data are consistent with other results published in literature; as an example, in fig. IS] we 
show the results that we obtained on a 8^ lattice at [3 = 0.96, which can be directly compared 
with the ones displayed in the top right corner plot of fig. 2 in [8], that were obtained for the 
same values of the parameters, but from simulations in the direct model. 

The data obtained from simulations on a 16^ lattice show no relevant discrepancy with respect 
to the 8^ lattice results for the values of (3 that we studied, except for the last one (/3 = 1.01), 
which is very close to ^c- this comes as no surprise, since, in that case, the typical length 
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0.14 



Figure 9: Comparison of results for the longitudinal component of the electric held obtained from 
simulations on a 8^ and on a lattice, at (5 = 1.01. For the 8“^ lattice, results are shown up to 
p/a = 4 only. 

scale given by the square root of the inverse string tension is of the order of some lattice units 
(roughly: ~ 4a), and the 8^^ lattice is expected to be affected by some relevant hnite size 

effects. Comparison of results obtained from the 8“^ and 16^ lattices at (5 = 1.01 is shown in hg. 

Testing the prediction of eq. (gH) is not a trivial task, since it is expected to describe the 
asymptotic behaviour of Ex at large distance from the interquark axis only, where the signal is 
decaying exponentially fast with /?; in order to explore large enough distances, we focused our 
attention onto the data obtained from lattices of size 16^. The results obtained from simulations 
at different values of coupling (ranging from /3 = 0.96 to /3 = 1.00) have a very good scaling 
behaviour, once they are expressed in terms of Sommer’s scale ro, as hg. shows: in this 
semi-logarithmic plot, the data appear to collapse onto a straight line, in agreement with the 
exponentially decaying asymptotic behaviour predicted by eq. (P|) . As a side remark, please note 
that fig. El also shows that the data enjoy good rotational invariance, since results corresponding 
to various off-axis directions display a behaviour which is fully compatible (within the errorbars) 
with those measured on the y or z axis. For the values of (5 that are considered here, tq roughly 
ranges from 1.8a (corresponding to (3 = 0.96) to 2.9a (for /3 = 1.00, in agreement with [10]). These 
results have been used to extract the dual gauge boson mass from the fit to the behaviour 
predicted by eq. gH); since the latter is expected to describe the asymptotic, large-distance, 
behaviour only, we discarded the data corresponding to distances shorter than a minimum value 
Pmin; tab-Elshows the results of the fit, for different values of Data analysis gives hints that 
a reliable choice for the lower distance cutoff to extract the dual gauge boson mass must be at 
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Figure 10: Behaviour of the longitudinal component of the electric held (displayed in a logarithmic 
scale) as a function of the distance from the interquark axis, for different values of /3 in the conhned 
regime; the results are scaled according to the corresponding values of Sommer’s scale tq (which 
was obtained from independent measurements). The dashed line is the result of the ht to eq. (I44|l . 
discarding the data which correspond to p/ro < 0.8. 
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Pmin/rO 

number of d.o.f. 

m^ro 

1.1 

11 

3.47(93) 

1.0 

14 

4.07(57) 

0.9 

17 

3.16(33) 

0.8 

18 

3.38(25) 

0.7 

21 

3.30(15) 

0.6 

23 

3.26(11) 

0.5 

26 

3.20(12) 

0.4 

29 

3.24(9) 


Table 6: Results for the dual gauge boson mass from the fit of eq. (EH) to the numerical 
data, for different values of the lower cutoff on distances 


least /Omin — 0.8ro: the corresponding curve is also shown in fier. llOl The results obtained with this 
choice appear to be in quantitative agreement (within two errorbars) with the ones published in 
the literature. 

A comment is in order: as opposite to the measurements of the inter quark force at large 
distances, in the study of the electric field profile our numerical simulations of the dual model 
did not account for particular advantages with respect to state-of-the-art numerical studies in 
the standard formulation of compact QED: thus in this case it seems that to improve the result 
precision it is necessary to resort to larger statistics. 

6 Conclusions 

In this work, we addressed a numerical study of compact U(l) LGT in L) = 4, a theory which 
possesses a confining phase analogous to non-abelian models. We restricted our attention to 
the infinitely heavy quark limit, and studied the behaviour of a confined QQ system, focusing 
our attention onto the interquark potential and force, and on the longitudinal component of the 
electric field induced in the mid-plane between the static charges. 

We discussed the analytical transformation mapping the original system to a dual model 
formulated in terms of integer-valued variables, which opens the possibility to perform efficient 
simulations. In particular, this method allows to overcome the exponential growth of the relative 
error at large interquark distances (which affects measurements of the interquark potential and 
force in the standard approach to LGT’s), focusing onto ratios of Polyakov loop correlation 
functions. To our knowledge, this is the first time this technique has been used to study these 
observables in compact U(l) LGT in D = 4. 

The numerical results that we obtained confirm the exponential error reduction with respect 
to Monte Garlo simulations in the conventional setting, and thus the method can be considered 
as a possible alternative to state-of-the-art error reduction algorithms, at least for the class of 
models for which this formulation is feasible. We measured data in a range of parameters which 
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includes and generalises the results obtained with modern techniques, finding agreement with the 
data available in the literature. Then we compared the numerical results with the theoretical 
predictions. For the interquark potential and force, we compared our data with the expectations 
derived assuming that an effective string picture holds in the IR regime: we summarised the 
assumptions underlying this scenario, discussing both the theoretical issues (including the aspects 
that have been pointed out most recently) and the problems that are still open in the comparison 
with numerical results. Our findings are in agreement with the conclusions of analogous recent 
studies, confirming the validity of the effective string picture at large distances. We also extended 
this study to the finite temperature setting, presenting quantitative results for the string tension 
in a range of temperatures up to about 3/4 of the critical deconfinement temperature. 

Finally, we addressed the study of the longitudinal component of the electric field induced 
in the mid-plane between the static sources, comparing our numerical data with the prediction 
from the dual superconductor scenario, and measuring the dual gauge boson mass which, in the 
confined phase of the theory, drives the exponential-like decay of the electric field at large distance 
from the interquark axis. 
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